Draft version February 2, 2008 

Preprint typeset using I£T^X style cmulatcapj v. 9/08/03 



CO 
O 
O 

CN 

> 

O 

O 

a\ : 

CN : 
> 

o\ : 

CN . 

■ 

ON 1 

o : 

cn . 
o ■ 

^: 

Or 

6 ■ 

u: 

&i . 



X 



SUNYAEV-ZEL'DOVICH EFFECTS FROM QUASARS IN GALAXIES AND GROUPS 
A. Lapi 1 , A. Cavaliere 1 , and G. De Zotti 2 

Draft version February 2, 2008 

ABSTRACT 

The energy fed by active galactic nuclei to the surrounding diffuse baryons changes their amount, 
temperature, and distribution; so in groups and in member galaxies it affects the X-ray luminosity 
and also the Sunyaev-Zel'dovich effect. Here we compute how the latter is enhanced by the transient 
blastwave driven by an active quasar, and is depressed when the equilibrium is recovered with a 
depleted density. We constrain such depressions and enhancements with the masses of relic black 
holes in galaxies and the X-ray luminosities in groups. We discuss how all these linked observables 
can tell the quasar contribution to the thermal history of the baryons pervading galaxies and groups. 
Subject headings: cosmic microwave background - galaxies: clusters: general - quasars: general 



1. INTRODUCTION 

The thermal energy content of the hot interstellar 
(ISM) or the intracluster medium (ICM) pervading 
galaxies or their groups and clusters can be probed with 
the Sunyaev-Zel'dovich (1980, SZ) effect. This arises 
when the hot electrons Compton upscatter some of the 
CMB photons crossing the structure; then the black body 
spectrum is tilted toward higher frequencies. 

In the [iwave band the tilt mimics a diminution 
AT^ W w — 5.5 y K of the CMB temperature, propor- 
tional to the Comptonization parameter y cx n T R. This 
is evenly contributed by the electron density n and the 
temperature T; in fact, what matters is the electron pres- 
sure p — nkT integrated along a line of sight (see Fig. 1): 



y = 2 



<7y 



d£p(r) . 



(1) 



To now, SZ signals have been measured in many rich 
clusters at levels y w 10 -4 or AT^ W sa —0.5 mK (see 
Rephaeli 1995; Birkinshaw 1999; Zhang & Wu 2000; 
Reese et al. 2002). These levels are consistent with ICM 
temperatures kT w 5 keV, sizes R of a few Mpcs, and 
central densities n ~ 10~ 3 cm" 3 . Similar values are indi- 
cated by the standard cluster view based on gravitational 
potential wells of virial depth kT v cx G M/R dominated 
by the mass M ~ 10 15 M Q of the dark matter (DM). In 
such wells ICM masses m w 0.2 M are settled in hydro- 
static equilibrium at specific energies kT ~ kT v . 

The above values also fit in with the X-ray luminosi- 
ties L x oc n 2 R 3 \[T w 10 44 - 10 45 ergs s" 1 emitted 
through thermal bremsstrahlung by the ICM. Groups, 
on the other hand, are underluminous relative to clus- 
ters; they emit far less than the baseline level L g cx T 2 
scaled at constant m sw 0.2 M after the DM rules. The 
observed Lx — T v correlation is clearly steeper, and goes 
from Lx cx T 3 in clusters to Lx oc T 4 or T^j in poor 
groups, albeit with a wide variance (O 'Sullivan, Ponman 
& Collins 2003). In other words, the ICM in groups ap- 
pears to be underdense compared to clusters. 
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The origin of such lower densities is currently debated. 
One view centers on extensive radiative cooling (Bryan 
2000) which would remove much low entropy gas. An 
alternative line of explanations (see Cavaliere, Lapi & 
Menci 2002, CLM02; and refs. therein) focuses on the 
energy injections affecting the ICM equilibrium while the 
DM is hierarchically accrued over dynamical timescales 
td- The inputs are provided when the baryons in galaxies 
condense to form stars possibly in starbursts, which then 
explode as SNe; alternatively or correlatedly (Menci et 
al. 2003, in prep.), the baryons accrete onto supermas- 
sive black holes (BHs) energizing active galactic nuclei 
(AGNs). Such feedback actions deplete the ICM density 
in the shallower potential wells by causing from inside 
thermal outflow and dynamical blowout; they also pre- 
heat the gas exterior to the newly forming structures, 
and so hinder its inflow. 

In any case, for groups in equilibrium the values of y 
can be anticipated from the continuum Lx through the 
model-independent relation (Cavaliere & Menci 2001) 



y/y g = {L x /L g ) 1/2 (T/T v f 4 



(2) 



Here y g cx (1 + z) 3 / 2 Tv is the baseline value scaled to 
the formation redshift z (Cole & Kaiser 1988). So for 
groups where Lx < L g holds we expect depressed y. 

Are enhanced SZ effects also possible, or even likely? 
What can these tell about the processes affecting n and T 
in groups and galaxies? Here we propose that a specific 
answer will come from SZ observations. 

2. THE TRANSIENT REGIME 

We start from recasting y g cx E/R 2 in terms of the 
gas thermal energy E cx pR 3 at equilibrium. A small 
group or an early massive galaxy with their virial tem- 
peratures kT v k 1 or 0.5 keV would produce SZ signals 
AT KU) /0.5mK w -5 or -3 x 10~ 2 (1 + zfl 2 . Larger en- 
ergies AE > E added to the ICM/ISM are expected to 
enhance the SZ signals yielding y/y g ~ AE/E. 

Such may be the case with the AGNs (see Valageas 
& Silk 1999; Wu, Fabian & Nulsen 2000; Roychowdhury 
& Nath 2002) that produce large total outputs, typically 
around 5 x 10 61 ergs over times around 10 8 yr comparable 
to td of their host structures. Such outputs can drive 
a blastwave sweeping out the gas mass and raising its 
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Fig. 1. — The geometry underlying Eq. (1). For a point in the 
structure r is the radial coordinate, s is its projection on the plane 
of the sky, and I is the coordinate along the line of sight. On 
the vertical axis we outline the initial density run, and the flow 
perturbed by the AGN-driven blastwave. 



pressure (see Platania et al. 2002; also Yamada & Fujita 
2001). 

We will see that the blast is constrained by the cou- 
pling of the AGN outputs to the surrounding gas, and re- 
strained by the initial pressure and the DM gravity. We 
describe the blast flow on using the self-similar hydro- 
dynamical solutions presented by CLM02. The simplest 
one obtains when the energy AE(t) oc t is delivered over 
times of order td at the center of an isothermal configura- 
tion for the DM and for the gas with n(r) oc p(r) oc r~ 2 . 
Then the leading shock moves out with uniform Mach 
number A4, i.e., with radius R s = R s t; the kinetic, the 
thermal and the gravitational energies of the perturbed 
gas all scale like R s . So we can consistently define inside 
R s the total initial energy E (modulus); most impor- 
tant, our solutions provide realistic predictions not only 
for strong blasts but also for the weak ones driven by 
constrained values of AE/E. 

Detailed profiles are presented in Fig. 2. Note that 
the perturbed flow is confined to a shell; this is bounded 
by the leading shock at R s , and by a trailing contact 
surface ( "piston" ) located at A R s < R s where the den- 
sity diverges weakly while the pressure is finite. So the 
relevant quantities may be also obtained from the simple 
and precise "shell approximation" (see Ostriker & McKee 
1988), which provides the energy balance in the form 

. „ 1 2 3 . , T , GMm 
AE-E=-mv 2 + -(p)V-^— . (3) 

Here (p) = p 2 (5 M 2 + 7)/ (5M 2 — 1) is the mean pressure 
within the shell volume V = Air R? s (1 — A 3 )/3, important 
for the SZ signals; M and m are the DM and the gas 
masses within R s . The Rankine-Hugoniot jump condi- 
tions yield the postshock quantities: the velocity v 2 = 
3R S (M 2 - 1) /AM 2 , the pressure p 2 = pi (5M 2 - l)/4 
and the density n 2 = AmM 2 /(M 2 + 3), given the 
preshock values pi , n\. 

As to the energy AE actually injected over td, the 
paradigm of supermassive BHs for the AGNs implies 
AE ps 2 x 10 62 / {Mbh/W 9 M q ) (1 + z)~ 3 / 2 ergs when 
the mass Mbh is accreted with conversion efficiency of 
order 10 _1 . The fractional energy / coupled to the sur- 
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Fig. 2. — Radial distribution of density and pressure in the blast- 
wave, relative to the post-shock values ri2 and P2 ■ Solid lines refer 
to a strong shock with AE/E = 3, dashed lines to an interme- 
diate shock with AE/E = 1, and dotted lines to a weak shock 
with AE/E = 0.3; the corresponding piston positions are given by 
A ss 0.7, 0.6, 0.45. 



rounding gas is poorly known. Including inefficiencies 
due to low momentum transfer, non-spherical geometries 
and covering factors it may range from / ps 10 -2 for 
radio-quiet, to some 10 _1 for strongly absorbed (BAL) 
or radio- loud quasars, a small minority. Average values 
/ ss 5 x 10~ 2 are shown below to be consistent with the 
observations both of the relic BHs in galaxies and of Lx 
in groups. 

The ratio AE/E is seen from Eq. (3) to be uniform, 
and to constitute a key parameter for the shock strength. 
Table 1 presents quantities relevant to our computations 
of SZ signals. 

We have extended our basic solution to initial density 
runs n oc r~ w , with 2 < u> < 2.5. Then the initial temper- 
ature declines as T(r) oc r 2 ~ u , and the shock decelerates 
as R s oc t 2 ^ , again under self-similarity; this also implies 
declining source luminosities with AE(t) oc £ 2 (5-2w)/w_ 
Table 1 shows that larger u> yields somewhat higher M. 
and (p) at given AE/E. 

3. TRANSIENT AND EQUILIBRIUM SZ EFFECTS 

In computing how y is enhanced during the blast tran- 
sit, we focus on y oc R~ 2 f dssy(s) averaged over the 
structure area, which will subtend small angles R/ Da ^ 
1' for an early group or galaxy. Normalizing the shock 
position as x = R s (t) / R, we find the full signal 



i = iii (1 _ A 3 )x+v /r— 

Vg 3pi 3pi V 



A 3 ), (4) 



in terms of y g = (4ctt /m e c 2 ) p(R) R. The last approx- 
imation applies for x ps 1, which maximizes the transit 
time in the structure and optimizes the observability. 

Strong SZ signals are seen from Eq. (4) and Table 1 
to require substantial blasts driven through the ISM 
or the ICM, i.e., input AE competing with the equi- 
librium value E. For w = 2 the latter writes E ps 
2 x 10 61 (ZcT/keV) 5 / 2 (1 + z)~ 3 / 2 ergs; so the ratio reads 



AE 
~E 



0.1 



/ M BH (kT v \ 
10- 2 1O 9 M \keV J 



-5/2 



(5) 
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Table 1. Relevant Blastwave Quantities 
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32.6 



Here Mbh is the mass accreted within td oc (1 + z)~ 3//2 
by the central BH in a massive galaxy, or by the sum of 
BHs shining within a group. Eq. (5) yields AE close to E 
for a poor group with kT v rts 1 keV and Mbh ~ 10 9 Mq. 
In going toward local, rich clusters Mbh clearly lags be- 
hind M, so AE/E will decrease strongly, see CLM02 for 
details. 

At the other end, toward galaxies AE/E is constrained 
not to exceed a few, lest the gas contained within kpes 
and the accretion it feeds are cut down (see Silk & Rees 
1998). The pivotal value AE/E w 1 recast in terms of 
the DM velocity dispersion a = (kT v /0.6 nip) 1 / 2 reads 

Converting to the bulge dispersion cr* oc a 1 ' 2 (see Fer- 
rarese 2002) yields M B h oc ct 4 . For values /«5x 10~ 2 
the relation accords with the observations in Tremaine 
et al. (2002). 

Our results are represented in Fig. 3 vs. the depth 
kT v of the host potential well. The square illustrates the 
minimal enhancement we expect from an early group 
at z = 1.5 with kT v = 1 keV, / = 5 x 10" 2 and 
Mbh = 1O 9 M , so with AE = 0.5 E. The bar gives 
a realistic upper bound for structures with steeper n(r), 
namely, with u> — 2.4; here E is larger but the en- 
ergy release is more impulsive, resulting (see Table 1) 
in stronger signals. With radii R « 250 kpc, the an- 
gular sizes 2R/Da ~ 1' are close to their minimum in 
the concordance cosmology (cf. Bennett et al. 2003). 
Comparable resolutions will soon be achieved, see § 4. 

The circles in Fig. 3 represent our results for a massive 
(a = 300 km s _1 , R w 100 kpc) and still gas-rich (to = 
0.15 M) protogalaxy at z = 2.5. The open circle refers to 
AE = E or Mbh = 6 x 10 8 M Q ; the filled one to AE = 
3E or M B h = 2 x 10 9 M o , just compatible with the 
scatter in the Mbh — o correlation. The related angular 
sizes are around 0.5'; with resolution fixed at 2 9b ~ l'j 
the signals will be diluted after (R/DaOi,) 2 ~ 1/4 and 
scaled down to AT^ W w — 20^iK. 

The inset represents the corresponding statistics. This 
is evaluated on inserting the related blue luminosities 
L = AE/I0ft d w 5 x 10 45 and 1.5 x 10 46 ergs s^ 1 (with 
a bolometric correction 10) in the quasar luminosity func- 
tion observed for z < 2.5 by Boyle et al. (2000), and dis- 
cussed by Cavaliere & Vittorini (2000, CV00). In terms 
of the cumulative fraction of bright galaxies hosting a 
type 1 quasar brighter than L, this reads 

N(L) L « 2 10- 2 (1 + zf' 2 (L b /L) 2 ' 2 , (7) 

beyond the break at L b = 5 x 10 45 [(1 + z)/3.5] 3 ergs s" 1 . 
The same luminosity function interpreted in terms of in- 
teractions of the host galaxy with its group companions 




kT v (keV) 

Fig. 3. — Predicted SZ signals as a function of the virial temper- 
atures of galaxies, groups and clusters. Dotted line: the baseline 
y g at 2 = 0. Shaded strip: central signals from gas in the equilib- 
rium set by SN preheating (Cavaliere & Menci 2001). Thick lines: 
same for feedback from AGNs with Mgj = 10 9 Mq and coupling 
/ = 3 X 10~ 2 (solid) or / = 10 — 1 (dashed). Square: area-averaged, 
undiluted signal from a group at z = 1.5, driven by AGN activity 
with Mbh = 10 9 Mq and / = 5 X 10~ 2 ; the bar represents the 
bound for u> = 2.4. Circles: same from a massive galaxy at z = 2.5, 
for AE = E (open) or 3 E (filled); the inset specifies the statistics. 



(CV00) yields a few signals per 10 poor groups, with the 
strength represented by the square in Fig. 3. 

After the passage of the blast, the gas recovers hydro- 
static equilibrium. This may be described by n(r) — 
n(R) exp(f3A<f>) for a nearly isothermal ICM in the 
(normalized) DM potential well A<fi(r), see Cavaliere & 
Fusco-Femiano (1976). The blast heats up the gas so de- 
creasing the value of the parameter (3 = T v /T, that we 
reset on using T averaged over the mass in the shell. The 
blast also ejects gas and depletes all densities; we reset 
n{R) by requiring the volume integral of n(r) to equal 
the gas mass left by the blast inside R at t — t^. 

The equilibrium SZ effect is then computed after 
Eq. (1), and the resulting signals are also plotted in 
Fig. 3. We recall from CLM02 that the equilibrium also 
provides good fits to the observed Lx in groups for the 
same coupling / w 5 x 10~ 2 , consistent with Eq. (6). 

4. DISCUSSION AND CONCLUSIONS 

This Letter centers on how the SZ effect is affected by 
the energy fed back by AGNs into the surrounding gas. 
We predict both transient, enhanced and long-term, de- 
pressed SZ signals to be produced by an energy addition 
AE to the equilibrium value E. 

In 1 keV groups the condition AE < E holds, and 
the AGN feedback has a considerable impact. A blast 
is driven through the gas; during its transit the area- 
averaged SZ signal is enhanced as the gas is just redis- 
tributed while its pressure is raised. But eventually a 
considerable gas fraction is ejected; so n is depleted and 
y cx n T is decreased at equilibrium. We have computed 
these effects under the restraints set to (weak) blast prop- 
agation by initial pressure and DM gravity. 

The result on the SZ effect is twofold: on scales 1' we 
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predict transient enhancements up to AT^ « — 80 /iK 
(a representative example at z « 1.5 is given in Fig. 3), 
followed by long-term depressions. The latter correlate 
after Eq. (2) with the equilibrium X-ray luminosities 
Lx oc n 2 y/T which are very depressed. 

Larger AE/E, yet constrained by Eq. (6), yield 
stronger SZ enhancements in gas-rich massive proto- 
galaxies with halo radii R ~ 100 kpc; representative ex- 
amples are shown by the circles in Fig. 3. With angular 
sizes 0.5', these may be diluted down to AT^ W as — 20 /iK 
when observed at a resolution around 1'. 

Such resolutions will be achieved by several instru- 
ments now being built or designed, enabling "blind" sky 
surveys for SZ signals to /iK sensitivities over tens of 
square degrees (see Carlstrom, Holder & Reese 2002). In 
particular, promising perspectives are offered by multi- 
beam, high frequency radio receivers like OCRA (Browne 
et al. 2000), and also by interferometers equipped 
with wide-band correlators like ATCA, SZA (Holder et 
al. 2000), AMI (Jones 2002), and AMiBA (Lo 2002). 
The SZ signals we consider may contribute equally or 
more than clusters to the excess power already detected 
at high multipoles with BIMA (Dawson et al. 2002). 
In the (sub)millimetric band the SZ signal is positive, 
and will be accessible to large bolometer arrays like 
BOLOCAM, whose developments will enable deep, wide 
surveys (Mauskopf et al. 2002). Eventually, ALMA 
(http://www.alma.nrao.edu/) will provide in selected 
areas higher resolution for both sides of the SZ effect. 

Enhanced signals as discussed here would constitute 
signatures of strong feedback caught in the act. This 
is specific of AGNs, since SNe feed back at most 0.3 
keV/particle (see CLM02); on the other hand, extended 
cooling which depletes n without increasing T hardly 
could enhance y oc n T. Interlopers might be introduced 
by merging events; however, these primarily govern the 
growth of the DM halos and set the virial T v oc M 2 / 3 
included in our baseline y g oc T^ 2 . Only an exceptional 
major merging may contribute an energy step sizeable 



but still bound by AE < E. Even this produces tran- 
sonic inflows in the high-T„ partner gas, originating lim- 
ited warmer features as picked up by highly resolved X- 
ray studies of clusters. Still smoother inflows are pro- 
duced by SN preheating (see Voit et al. 2003), while 
stronger blasts are driven by AGNs, in the galaxies and 
groups that we propose here as primarily SZ objects. 

Detecting 10 such signals will require surveys over 500 
arcmin 2 at 1' resolution, based on the conservative sur- 
face density of 10 2 powerful quasars/deg 2 consistent with 
Eq. (7). For groups our evaluations (Fig. 3) lead to SZ 
enhancements in the range y/y g ~ 1.2 — 4 over a de- 
pressed if wrinkled landscape. In fact, for kT v « 1 keV 
the baseline y g is affected mainly by SNe (see Fig. 3); 
these depress the average levels somewhat, and cause at 
z « 1.5 a 10% relative scatter (at 96% confidence). This 
landscape may be sampled or bounded from independent 
groups catalogued at comparable z. For massive proto- 
galaxies intrinsically stronger enhancements, less depres- 
sion and narrower scatter obtain. The candidate peaks 
are to be followed up with ALMA for higher resolutions; 
in addition, optical z, and optical a or X-ray T will re- 
quire current or moderately extrapolated techniques (cf. 
Rosati, Borgani & Norman 2002; Shields et al. 2003). 

In conclusion, we expect that AGN energy outputs 
around 10 62 ergs with coupling / ~ 5 x 10~ 2 leave two 
consistent relics: the depressed X-ray luminosities Lx 
in local galaxies and groups (see CLM02); the Mbh — o~ 
relation on subgalactic scales (Eq. 6). Relatedly, on in- 
termediate scales we evaluate here (Fig. 3) transient SZ 
signals standing out of a generally depressed landscape. 
Such signals can provide real time evidence of AGN feed- 
back acting on the diffuse baryons in galaxies and groups. 
The evidence should be looked for primarily in the SZ 
surveys that will be soon available. 



We thank N. Menci and our referee for helpful com- 
ments. 
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